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I. INTRODUCTION 

The measurement problem in quantum mechanics remains, almost a century after the theory's formulation, the 
major obstacle to considering the theory as truly fundamental. Despite heroic efforts by many insightful physicists, 
the difficulties involved have not yielded and wc are still lacking any fully satisfactory option. The basic issue, as 
described for instance in Q, is the fact that the theory relies on two different and incompatible evolution processes. 
Using Penrose's characterization^, there is the U (unitary) process, where the state changes smoothly according to 
Schroedingcr's deterministic differential equation, and the R (reduction) process, in which the state of the system 
changes instantaneously, in an indeterministic fashion. The U process is supposed to control a system's dynamics all 
the time that the system is left alone, while the R process is called upon whenever a measurement has been carried 
out. 

The problem is that no one has been able to characterize in general when a physical process should be considered 
a measurement. This issue has been studied and debated extensibly in the scientific and philosophical literature [f|. 
Of course, in laboratory situations, one clearly knows when a measurement has been carried out. Nonetheless, as 
characterized by J. Bell @, a for all practical purposes (FAPP) approach is not satisfactory at the foundational level, 
as it involves treating the system differently from the measuring device or the observer, and this division is one for 
which the theory offers no specific internal rules. Its most conspicuous inadequacy occurs in cosmological applications, 
where one can not use any interpretation that relies on an observer, or on a measurement device. [7[, 

One approach^ to resolving this problem of standard quantum theory is to modify it by incorporating novel 
dynamical features responsible for 'the collapse of the wave function'. It may be characterized as the promotion 
of quantum theory from a theory of measurement to a theory of reality, in which some of the physical properties 
of a given system take values, regardless of whether they are observed or not. Such an approach can be applied 
to cosmology, and we shall do so by focusing on the problem of emergence of the seeds of structure in inflationary 
cosmology. 

According to the standard account of inflation, in its early stages, a relatively generic state of the univers^H is 
driven towards a homogenous and isotropic Robertson Walker space-time. This expands almost exponentially, driven 
by the potential of the inflaton field, which acts as a large effective cosmological constant. The inflaton field is a 
quantum ficlc0 which is taken to initially be in the so called Bunch-Davies vacuum state. This state is completely 
homogeneous and isotropic, and the dynamics preserves these properties. There are various proposals to address this 
issue. One is to suppose that the state vector can be expanded in some 'natural' basis and, somehow, it is one of 
these basis states which describes our inhomogeneous and anisotropic universe. As discussed in detail in [J], none of 
them can be regarded as satisfactory. 

There are analogous instances where quantum theory presents us with a symmetric quantum state whereas nature 
exhibits asymmetric behavior. Perhaps the best known example is the a decay of a J = nucleus. Although the wave 
function is rotationally invariant, the alpha particles are seen to move on linear trajectories. The problem was studied 
by N. MottfH, and was resolved by heavy usage of the collapse postulate appropriate to a measurement situation[l|. 
However, in the cosmological problem at hand, even if we wanted to, wc could not achieve a similar explanation, for 
we can not call upon any external entity making a measurement. 

The quantum fluctuations of the inflaton field, according to the standard accounts (which are phenomenologically 
quite successful) are supposed to seed the growth of structure and lead to the formation of galaxy clusters, galaxies, 
stars, etc. There is no way we can invoke anything like a measurement postulate as part of the explanation of the 
emergence of those primordial seeds of structure. Not only there were no observers or anything to play the role of 
a measuring apparatus at the time but, because we and our measuring devices (and indeed any conceivable kind of 
astronomers of alien civilizations) owe their existence to the process generating those seeds, they can not be any part 
of their cause! 

To summarize, the problem is that one cannot explain the emergence of the observed asymmetries in homogeneity 
and isotropy when one has a theory with an initial state, the Bunch-Davies vacuum on the RW space-time, which 
is 100% homogeneous and isotropic and a dynamics that does not break such symmetries. Nonetheless, we wish to 
recover the phenomenological success of the standard account, where quantum fluctuations in the vacuum state are 
the source of the first inhomogeneities and anisotropics 

In the case of alpha-decay, dynamical collapse theories, either the GRWfIc| model or the CSL0 model, satisfactorily 
resolves the problem. In CSL, which is usually considered to have superseded GRW, the Schrodingcr equation is 



There are several conditions that are required for this, but we shall not elaborate on that here. 

2 There are scalar, vector and tensor fields, but we will focus here on just the scalar field which is responsible for the anisotropics that 
have been observed in the CMB, so far. 

3 The fact that, due to the kinetic term in the energy momentum tensor of the scalar field, the space-time can not exactly correspond to 
the de Sitter metric, implies that the state is not exactly the Bunch-Davies state. However, this difference is negligible for our treatment. 

4 We emphasize that the so-called quantum fluctuations are quantum uncertainties and should not be confused with thermal or statistical 
fluctuations. 
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modified by adding to the usual hamiltonian a non- unitary term iHc, where Hq is a unitary 'collapse hamiltonian.' 
He depends upon a randomly fluctuating classical field w(x, t). The eigenstates of He are essentially mass density 
cigenstatcs. This modified Schrodinger equation is supplemented by a second equation, the Probability Rule, which 
gives the probability that nature chooses a particular w(x, t). The dynamics is such that, for each u>(x, t) of high 
probability, the collapse hamiltonian evolves the state vector toward one or another eigenstate of mass density, and 
this occurs with the Born probability. 

In the description of alpha decay, under the usual hamiltonian, the state vector describes the alpha particle inter- 
acting with the gas atoms in its path. With the added CSL dynamics, as more and more gas atoms get involved, 
for each high probability u>(x, t), the collapse hamiltonian more and more rapidly drives the state vector toward an 
approximate mass-density eigenstate describing the alpha moving in a straight-line path among its associated atoms. 

Returning to the standard inflationary scenario, since there is no physical mechanism that could account for the 
asymmetries, it is natural to consider addition of a physical mechanism that can do so. The mechanism of wave function 
collapse in the cosmological setting was suggested and modeled in an ad hoc and phenomenological way, in previous 
work along this line.[ll[ This led to phenomenological constraints on the parameters characterizing the suggested 
collapse mechanism, in order to ensure the theory is able to reproduce the observational results. In this manuscript, 
we consider this idea anew, employing however the rather well developed CSL formalism for the description of wave 
function collapse. 

The article is organized as follows. Sections IHUfVl present the non-collapse setting, and the remaining sections 
invoke the collapse. 

In Section UH we discuss the inflationary paradigm, and give estimates of important numerical quantities. Section 
IIIII obtains the hamiltonian for the quantum perturbations on the inflaton field and presents the dynamical solu- 
tions without collapse. Section IIVI discusses the the metric perturbation on the Robertson Walker metric known as 
the Newtonian Potential. The connection between classical gravity and quantum variables is obtained by invoking 
semi- classical gravity, whereby the Newtonian Potential is related to a certain quantum expectation value, and this 
important assumption is discussed. Section IVl reviews the relation of experimentally observed quantities to the New- 
tonian Potential, and thereby to those quantum expectation values. Section I VII then reviews the CSL formalism. 
Section IVlII adds the CSL modification to the hamiltonian evolution of the inflaton perturbations. Sections IVIIII and 
IIXI show how it is possible to obtain agreement with the observations. Section [X] uses this result to obtain expressions 
for some physical quantities, and their probabilities of realization. Sections IXII and IXHI discuss our conclusions. 

Regarding notation, we will use signature ( — h H — h) for the metric and Wald's convention for the Riemann tensor. 

II. INFLATION 

The starting point for the discussion of inflation is the action of a scalar field coupled to gravity. 

R\g] - l^a^ b <Pg ab - v[4>] (l) 

The field equations derived from the action Eq.([T]) are, 

G b l = 8nGT b a (2) 

where T b a is give by: 

T b a = g ac d c <pd b <j> + 5 a b {- l -g cd d c <pd d <P - V (</>)) (3) 

Now we address a fundamental issue, how to combine quantum theory and gravitational physics. Despite important 
advances, it is well known that we are still lacking a fully workable and satisfactory theory of quantum gravity. It 
is also well known that a quantum gravity approach based on canonical quantization leads to what is, in effect, a 
timeless theory. [l2[ 

On the other hand, time seems not only an important aspect of any discussion of cosmology but, also, time is 
needed for the CSL dynamical reduction theory which we wish to incorporate. These considerations lead us to use 
the approach based on scmiclassical gravity, where matter fields are treated at the quantum level whereas gravitation 
(although quantum at the fundamental level) is assumed to be in a regime where it can be treated in a classical 
manner to a very good approximation. [l| The fact that inflation is thought to occur at energy scales well below the 
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Planck mass lends support to this assumption. (For further discussion on the issue and on the kind of treatment 
capable of incorporating dynamical reduction in such a context, see [l3lph This approach differs from that followed in 
standard treatments. There (as here), the starting point is a classical background for both the gravitation and matter 
fields, but a quantum treatment is used for the perturbations of both the metric and the inflaton. In the present 
work we adopt in principle the strict semiclassical approach described above, limiting the quantum treatment to the 
perturbation of the inflaton field. This difference should be kept in mind when comparing the standard approach to 
our treatment using the CSL theory. 

We now proceed to describe the basic setting of the problem, starting with characterization of the background 
metric and inflaton field, followed by treatment of the corresponding perturbed quantities. 



A. The Background. 

The analysis here will be based, as is usual, on separating the metric and scalar field into a spatially homogeneous- 
isotropic background part and a fluctuation part. That is, the scalar field is written <f> = 4>q + 5<f>, while the metric is 
written as g = go + Sg, where <po is a function of the conformal tim(J^| r\ only, and go characterizes the spatially flat 
Robertson Walker cosmology, i.e., 

dsl = a 2 (ri)[-dr] 2 + 8 lJ dx l dx : >}, <p (rj). (4) 

The background metric and scalar field are treated classically. The scalar field and Friedman evolution equations 
are 

/ o ' + 2^ / o + a 2 ^ o y[0o]=O, 3^ =4^G[^ 2 + 2o 2 V[^o]], (5) 

where the prime denotes derivative with respect to the conformal time rj, whose range during the inflationary era is 
negative, i.e., rj <G (—00, —6) with 6 > 0. 

The scale factor corresponding to the inflationary era of standard inflationary cosmology which follows from Eqs.© 
is, to a good approximation, 

a(ri) = (j^) 1+e - (6) 



where the Hubble constant (the expansion rate in time t) is 

Hi 



(87T/3)GV[cbo} m ■ (7) 



Written in terms of H = a' (rj) / a(rf) , the so-called slow-roll parameter e = 1 — %' /Ti 2 is considered to be very small: 
e<l during the inflationary stage. As for the scalar field </>o, the slow roll regime corresponds to 4>' = —{a 3 /3a')d<pV^\ 
Also, it is assumed that the change in the scalar potential is so small that the Hubble parameter Hj is essentially 
constant. 

According to the standard scenario, this inflationary era is followed by a brief reheating period in which the 
universe is repopulatcd with ordinary matter fields (a process that, for simplicity, we take to occur instantaneously), 
immediately followed by the radiation-dominated era, which begins the evolution of standard hot big bang cosmology, 
leading up to the present cosmological time. While the functional form of a(rj) changes during these later periods, 
that is irrelevant for our calculation, which deals solely with the inflationary era where it is assumed that the seeds 
of structure are generated. 

We shall take inflation to start at rj — — T, assume that the inflationary regime ends at the start of the radiation- 
dominated era at r\ = — t, and set a = 1 at the 'present cosmological time.' 

The effects of the late time physics, comprising the physical processes occurring between the time of reheating and 
the time of decoupling (where the matter becomes transparent to the radiation), to the extent that we shall refer to 



In contrast with the standard approach where the quantum-classical cut is tied to the background-perturbation separation, in the 
treatment developed in reference|l3||. the quantum-classical cut is tied to the separation of gravitation and matter fields. In particular, 
the so called zero mode of the inflaton field is treated at the quantum level. 

Note that a(r))dr) = dt, where t is ordinary time, so if a(?y) ~ — then a(rj) grows exponentially with t. 

This slow roll stage, corresponds mathematically to the "terminal velocity" of a body subject to a constant force in addition to a friction 
term. Thus the condition is = which corresponds in conformal time to <f>Q = Hcpg. Using this in the first equation in J5} leads to 
4>' = -(a a /3a')d,t,V . 
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them, will be taken as fully codified in what are called transfer functions (which relate the primordial fluctuations 
to those fluctuations that are directly observable in the CMB) and which will be, for the most part, ignored in the 
present work. 

The only additional information (beyond the behavior of quantities during the inflationary regime) we need, in 
order to make the estimates below, is the relationship between the scale factor at the end of inflation, the scale factor 
at the time of decoupling and the scale factor today. 

Of course, the quantities directly observed are those that were present at the time of the decoupling (whose conformal 
time is tjd) which lies in the matter-dominated era. However, we shall evaluate such quantities at the end of inflation, 
at r) = —r. We can do this because the transfer functions allow us to go backward from r/ = r/jj to r/ = — r. For 
example, the famous acoustic peaks, the most noteworthy feature of the CMB power spectrum, which are related to 
aspects of plasma physics, can be effectively subtracted out using the transfer functions. The results lead back (as we 
shall discuss later in detail) to the scale-free Harrison- Zeldovich spectrum at the end of inflation. Therefore, it is the 
H-Z spectrum for which our calculation must account. Alternatively said, if the transfer functions were constants, we 
would be directly observing the H-Z spectrum today. 

B. Estimates. 

Now, we need to estimate the values of the conformal time rj = — r at the end of inflation and the conformal time 
rj = —1~ at which the inflation starts. This may be done as follows. 

Recall that the temperature of radiation, regardless of era, scaled like 1/a. The radiation temperature today 
corresponds to 2.7 °K = 2.4 x 10~ 13 GeV. We adopt the widespread assumption that the inflation scale, i.e., the 
effective temperature at the end of inflation, corresponds to the Grand Unification Theory (GUT) scale of about 
I0 15 GeU. Therefore, we estimate a(— r) = 2.4 x I0 -28 . Next, we must find the value — r corresponding to this scale 
factor. 

Using the GUT inflationary scale for V ^(GUT scale) 4 , and as we use c = fi = I , G = Mp, , we employ Eq.([7]) 
to determine Hi m 3(M% UT /M Pl ) = 3(M GUT /M P i) 2 x 10 19 GeV. = 3 x 10 n GeV. Knowing a(-r) and Hj, we find 
t ~ fO -22 Mpc from Eq.((l|. (The result can be written r ~ {Mpi/Mqut) 10 -26 Mpc, where Mpi is the Planck mass). 

Next, in order to estimate the value of the conformal time at the start of inflation, rj = —T, we assume that 
inflation lasts, say, 70 e-folds (usually considered a lower bound for inflation to solve the naturalness problems). Thus 
T/t = a(—T)/a(—T) = e 70 ps fO 30 . Combining this with the previous result r ~ fO -22 Mpc. gives T = fO 8 Mpc. 
Note that this time becomes larger if we require inflation to last more than that number of e-folds. 

Also, we shall need an estimate of the co-moving wave number k = |k| for the modes k which are relevant for the 
observed CMB data. 

These are obtained by noting that the surface of last scattering/decouping is at cils = 2.7A'°/3000if ° = f/f fOO and 
that the co- moving radius R of the last scattering sphere is determined by the requirement that the photons that started 
from a point on that sphere at the time of decoupling are just reaching us today. This gives R ps 2/Hq ps 6 x f 3 Mpc. 
The physical radius of that sphere is then Rph ys = o-lsR ~ 5.5Mpc. 

Now, the scales that are relevant range from those corresponding to this radius to, say, 10 -5 of this value. This 
corresponds to angular scales of 2ir x f CU 5 which is the smallest scale that can be seen with current technology in 
the CMB. Thus the relevant modes are those whose physical wavelength aLS^Tr/k at last scattering roughly ** lies 
between R Phys and 10~ 5 Rp hys . This gives lO^Mpc^ 1 < k < f0 2 Mpc _1 . 

We emphasize that this discussion about the values of k is completely independent of the precise functional form 
of a(r]) between the end of inflation and the decoupling time: the only relevant data is the relationship between the 
values of the scale factor at the end of inflation and the scale factor at the decoupling time. 

It is clear that the conditions kT >> 1 and kr << 1, which we shall use to make approximations, hold by very 
wide margins for the relevant modes. 

III. THE FLUCTUATION OF THE INFLATON FIELD. 

Now we consider the fluctuation of the inflaton field, 5(j>(x,rf). We start with the perturbed action up to second 
order in the scalar field fluctuation, written in term of the auxiliary field y = a8(p, 



For a brief period prior to the end of reheating (which we take here to occur essentially instantaneously at the end of inflation), this 
relationship between the temperature and the scale factor does not hold, because the radiation and particle content of the universe are 
replenished at that time as a result of the process of inflaton field decay. 
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5S^ = I J dv d 3 x (V 2 - (Vy) 2 + U 2 y 2 - 2Uyy^j . 

The Lagrangian density is then, 

SC^ = ^y' 2 -(Vy) 2 +n 2 y 2 -2nyy' 

The canonical momentum tt conjugated to y is n = d6£W/dy' = y' - Uy. Note that S<t>' = n/a. With H 
a 1 \rj) I a(r\) and a(rj) = —1/Hir), the non-vanishing equal-time poisson bracket and the hamiltonian are: 



(8) 



(9) 



[y(x), tt(x')] = 5(x - x'), 6H {2) = - I 



7r 2 (x)-- 7 r(x)y(x) + (V 2 /(x)) 2 



(10) 



where here, and in what follows, wc suppress the dependence of all variables on r/. 
We next focus on the individual modes of the field: 



y( x ) 



(2 



J dky(k)e^, tt(x) = ^-L^ J dkn(k) t 



(11) 



In terms of y(k), 7r(k), which are no longer real, the non-vanishing equal-time poisson bracket and the hamiltonian 
then become 



[y(k),7r(k')]=5(k-k'), 8H^ = \ J dk 



7r(k)7r*(k) - -[7r*(kWk) +7r(k)2/*(k)l + k 2 y(k)y*(k) 



(12) 



In making the transition from this classical description to the quantum description, we replace a c-number real 
variable a by a hermitian operator d, and replace the poisson bracket by i^ 1 x (commutator bracket). (We shall work 
in the Schrodinger picture, so the operators will be time- independent, and the state vector, evolving according to 
Schrodinger's equation, describes the time evolution). Thus, j/(x)t = y(x), and therefore y(k)^ = y(— k). Similarly, 
7r(x)^ = 7r(x), and therefore Tr(k)^ = tt(— k). Therefore, because the classical variables describing the modes are 
not real, the operators y(k) and 7r(k) are not hermitian. In order to work with hermitian operators, we turn to a 
description in terms of symmetric and antisymmetric fields. 



A. Symmetric and antisymmetric fields. 

Still in the classical theory, we write each field as the sum of symmetric and antisymmetric parts: 



y( x ) = 7}[y( x ) + y(- x )} + 7^[y( x ) ~ y{~ x )] = ys(*) + jm(x), 



and similarly for 7r(x). Putting (fT3|) into (fTUj) . and using J dxfs(x)gA(x) = 0, Eq. (fTUf becomes 



(13) 



SH^ = - I 



tt|(x) - -7r s (x)ys(x) + (Vys(x)) 2 







+1 2J 


f dx 







- -7Ta(x)j/a(x) + (V2m(x)) 2 



(14) 



i.e., 6H& = +8Hf. 

The equal-time poisson brackets of the symmetric and antisymmetric fields follow from (|10l) : 

[y 5 (x),^ s (x')] = ±[<5(x-x')+<$(x + x')], [y A (x),7r A (x')] = l[<J(x - x') - <S(x + x')], 



(15) 



and [2/5(x),tta(x')] = [?m(x), tt s (x')] = 0. 
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B. Modes, symmetric case. 

Now, we consider the individual modes of the field, 
2/s(x) 



(16) 



Because ys(x), 7rs(x) are symmetric, we get k) = ys(k), tts(~ k) = 7Ts(k). Because j/s(x),7Ts(x) are real, we get 
2/5 (k) = 2/s(— k), 7rJ(k) = 7T5(— k). So, j/ s (k) = j/s(k), 7r s (k) = 7r5(k): these classical variables are real. Therefore, 
the non- vanishing poisson bracket is, using <|15[) : 



{ys(k),TT S (k')} = {y s (k),n* s (k')} = 



(27T) 



x axe 



[ys(x),7rs(x')] 



= -[<5(k-k') + 5(k + k')]. 



(17) 



In evaluating the hamiltonian, the first term is 

i | dx7r|(x) = \ J dx7r s (x)7r5(x) = J dx J dk J *V k - x e- ik '- x 7T5(k)7r5(k') 

= i y rfk7r|(k) = | rfk7r|(k). 



(18) 



In the last step, the integral over all k is converted to the integral over the upper half k-plane, since the lower half 
k-planc makes the identical contribution. Similar steps are to be taken for the other terms in the hamiltonian, which 
then becomes: 



SH 



(2) 



dk 



7ri(k)--7r s (k)y s (k) + A;^(k) 



(19) 



Finally, in limiting to the upper half k-plane, we note that the <5(k + k') term in the poisson bracket isn't used in 
calculating the equations of motion, so the poisson bracket is effectively just the first term in (fl7]) . However, there is 
a factor of 1/2 in the poisson bracket which makes this not quite canonical. So, we shall define new variables 



X s (k) = V2y s (k), Ps(k) = \/27T5(k). 
In terms of these variables, the equal time Poisson bracket and hamiltonian are 



(2) 



[X s (k),P s (k'))=S(k-k'),SH^ 
Lastly, we proceed to quantize, so Eq. ([2"TT) becomes 



dk 



P£(k)--P s (k)X s (k) + k'X^) 



z<S(k - k'), 5Hf = 1 I (k) - -[P s (k)X s (k) + X t 



(20) 



(21) 



(22) 



Again, we emphasize that, while the classical variables are functions of conformal time 77, we are choosing to work in 
the Schrodinger picture. Thus, the operators are 77-independent (except for the Hamiltonian where the //-dependence 
is explicit): the 77 dependence is relegated to the behavior of the state vector. 



C. Modes, anti-symmetric case. 

Here we proceed in a manner exactly parallel to that used in the previous case. 
We start with, 

2M(x) = J dke**iy A (k), 7r A (x) = J dke^n A (k). (23) 

However, as ?m(x), ^(x) are anti-symmetric, we get y A {— k) = — y A (k), 7r^(— k) = — 7T^(k). Because y A (x), ir A (x) 
are real, we get y* A {k) = — y A (— k), 7r^(k) = — ir A (— k). So, y* A (k) = y A (k), 7r^(k) = ^(k): the i-factors in the 
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definitions of y A (k), ^(k) were chosen to make these real so that they become hcrmitian operators in the transition 
to quantum theory. So, although 7r(k) is not hermitian, it has been expressed in terms of hermitian operators: 



7r(k) = 7r s (k) + iir A (k). 
The non- vanishing equal time poisson bracket is, using (|15j) : 



(24) 



[y A {k),K A {k!)] = [y A (k)^ A (k!)} = 



(2tt) 3 



x /ax e 



e lkx [j/ A (x),7r A (x')] 



-[<5(k-k')-5(k + k')]. 



(25) 



The rest of the argument goes through just as for the symmetric case. The only difference is that the sign of <5(k+k') 
is positive for the symmetric case, Eq. (fTTj) . but negative for the anti-symmetric case, Eq. (|2"5)) . However, when we 
limit to the upper half k plane, <5(k + k') plays no role. Therefore, after quantization, we obtain the anti-symmetric 
version of Eq. 



[X A (k),P A (k')] =i(J(k-k'), 5H A 2) = ^ I dk P|(k) - -[P A (k)X A (k) + X s (k)P A (k)] + k 2 X 2 A {k) 



1 



1 



V 



(26) 



In this way, the field has become a simple collection of independent modificd-harmonic oscillators, with each mode 
evolving independently. If we define 

X = VdkX a (k), P = VdkP a (k), (27) 

with indices a = S, A, k suppressed, the commutation relations and hamiltonian for a mode characterized by k = \k\ 
is 



[X, P] = i, H k 



P 2 - - [PX + XP] + k 2 X 2 



(28) 



D. Non-Collapse Theory: Expectation Values 

We emphasize that, in our treatment, the hamiltonian Hk is not the sole cause of the dynamics. We have yet to 
incorporate the CSL modification of quantum theory (Section I VII and beyond). Nonetheless, it is interesting, and 
shall prove useful, to calculate expectation values due to the collapse-free dynamics characterized by H k alone. 

This is easy to do. For any operator A, (ijj, t\A\ip, t) satisfies 



— (ip,r)\A\iP,ri) = -i{i),v\[A,H k ]\^,r]). 
drj 



(29) 



Because H k is quadratic, the set of expectation values of any power of operators form a closed set which can be 
solved. The initial condition is that the mode wave function is in the Bunch Davies vacuum, which is just the 
harmonic oscillator ground state, at the initial time rj = —1~ 



W>-T) 



-v 2 ' 2 \ (a#,-T) = (Tr/fc) 1 /^- 2 ^. 



(30) 



(^fc) 1/4 

Writing (A) = (ip, r/lAlip, ?/), the first order equations, the consequent equations of motion and their solutions arc 

(31a) 
(31b) 
(31c) 



drj rj 



± { P)=-k 2 (X) + ^- 
drj rj 



(X) = C 1 —e 



V 

ikif 



drf 



I ' 






krj. 


k 


krj. 



(P) = de lk, i +C 2 e- lk11 . 
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(It is notable that (P) has the usual harmonic oscillator solution, even though the Hamiltonian is not the usual 
harmonic oscillator hamiltonian.) From the initial conditions Eq. ([30|) . we see that that (X) and (P) vanish initially, 
C\ = C2 = 0, so for all 77, 



(X) = (P) = 0. 

The second order equations, with Q = (X 2 ), R = (P 2 ), S = ([XP + PX]), are 

jLq = S -^-, ^R = -k 2 S + —, ^S = 2[R-k 2 Q}. 
drj 7] drj r\ drj 



(32) 



(33) 



Because the algebra of the commutator brackets is identical to that of the poisson brackets, and for the classical 
variables the product of two solutions is the solution for the product, the same is true for the solutions of Eqs. ([33l) . 
They are the product of the solutions Eqs. (|31cj ): 



krj 



2 1 
-C 2 — e' 2lk ^ 



R = Cie 2ikr > + C 2 e- 2ik71 + C 3 , 



S= -2id\e 2ikv 
k 



krj 



k 2 



1 



krj 



2iC 2 -e- 2ikri 



1 - — 



1 

krj 



2 1 
2 



1 



1 



(krjf 



a 



'k 2 rj' 



The initial conditions are 

Q(-T) = l/2fc, R(-T) = k/2, S(-T) = 0. 
Assuming kT << 1, we obtain C\ = —C2 = 0, C3 = k/2 and so 



Q = (x 2 ) = 



1 r 

2k 



(krj) 2 ' 



, R = (P 2 ) 



S= ([XP + PX]) 



1 

krj 



(34a) 
(34b) 
(34c) 

(35) 
(36) 



It shall be seen that the result (P 2 ) = k/2 which, together with Eqs.flM]), (HOJ), (HZl),(j32) et. scq. implies 

(^(k)^(k')*) = ((7T S (k)+i7r A (k))(7r S (k / ) -iTTA(k'))) = (*s(k)7T S (k')) + (7T A (k)^(k')> 

2 

)6 kk/ = (P 2 )6(k-k r ) = ^S(k-k'), 



'2dk 



(37) 



in the usual treatment is what is cited as causing agreement between the theory of inflaton perturbations and the 
effectively observed ** H-Z spectrum of temperature fluctuations. 



IV. THE FLUCTUATION OF THE METRIC: NEWTONIAN POTENTIAL. 



In order to connect with observations, we need the metric perturbation ^ known as the Newtonian potential at 
the end of the inflationary period. As we have mentioned, immediately following that time, the universe entered 
into the brief so-called " reheating period" , when the inflaton field is supposed to have decayed, in a not completely 
understood manner, into the matter content of the present universe: dark matter, baryons, electrons, photons, etc. 
As we shall see, the Newtonian potential is related to the fluctuation (inhomogeneity) in the inflaton field, so that it 
ends up being tied to the inhomogeneities of the matter created during reheating. During reheating, it is supposed 
that the Newtonian potential did not change in any appreciable manner from its value at the end of inflation. Thus, 
this primordial Newtonian potential can be use to determine the ensuing evolution of the matter inhomogeneities, 
entailing well-known physics such as baryon acoustic oscillations, and other processes described by the so-called 
transfer functions. 

Eventually, the universe evolved to what is variously called the time of last scattering or the time of decoupling 
of the photons from the plasma, or the time of recombination of electrons and protons, when the atoms formed and 
universe suddenly became transparent to the radiation which we now detect. At that time, the universe is considered 
to have been in local thermal equilibrium. The quantity that we presently measure is the temperature variation as a 
function of coordinates on the celestial sphere, AT(8, (p)/T ( AT(9, ip) = T(8, ip)—T, where T is the mean temperature 
over the sky) . This temperature variation is due to the inhomogeneities in the matter density distribution at the time 
of last scattering. 
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The satellites measure the temperature by detecting the blackbody microwave radiation at a number of frequencies. 
If v is the frequency at the peak of the spectrum, we have the relations 

AT 5V 1 T 

-jr = ~ « 3*. (38) 

where the last step uses the transfer functions to strip away the physics ensuing between the end of reheating and the 
time of decoupling, effectively considering that we are directly observing ^. The effect of Vl 7 is, first, a gravitational 
red- or blue-shift. Second there is an effect on the rate of expansion of the universe whose combined (Sachs- Wolfe) 
effect gives the 1/3 factor @ . 

The Newtonian potential is related to the fluctuation of the inflaton field as follows. When the perturbation of the 
Robertson Walker space-time is taken into account, with the appropriate choice of gauge (conformal Newton gauge), 
and ignoring the vector and tensor part of the metric perturbations, Eq.(|3]) is replaced by 

ds 2 = a(r]) 2 [-(l + 2<^)dr, 2 + (1 - 2$>)6 ij dx i dx j ], 

where $ and \E' are functions of the space-time coordinates 77, x l . 

Next, consider Einstein's equations to first order in the perturbations. The expression for the energy-momentum 
tensor T b a for the inflaton field is Eq.Q. Its linear perturbation components are: 

5T ° = a- 2 [<p'<?$ - <f>' 8<t>' - d$Va 2 d(j)}, 5T? = d^-ar 2 ^^), 8T] = a _2 [#,<ty' - <t>'^ - d+VaHftf). (39) 

Then, Einstein's equations to first order, SG = 8irG5T, lead to = $ and 

V 2 * + /i* = AttG(lo54> + <jf Q 5<t>') (40) 

where fi = H 2 -W and u = 3H(p' + a 2 d^V . 

As discussed following Eq.(|6]), the slow-roll approximation corresponds to lo = 0. We are ignoring here terms of 
order e, which implies /i = 0. Thus Eq. (|40[) and its Fourier transform become: 

V 2 *(7 ? , x) = 4^^(77)^' (77, x) = 47rG? y ^ 77(77, x), -fc 2 *(77, k) = 4^^(77)^(77, k) = ^^K fa k). (41) 

Now we make the transition to quantum theory. The usual procedure is to quantize both sides of Eq. (|41[) . so that 
#(k) — 7r(k). But then, as a consequence of Eq.flU]) and Eq.([32]). (*(k)) ~ (^(k)) = 0. Therefore, (ip,r}\if(x)\il),Ti) = 
0. 

How is this to be interpreted? One would like to identify the expectation value of the Newtonian potential operator 
with the value of the Newtonian potential in nature. However, for a state \ip,v) representing nature, the expectation 
value should vary with position. Therefore, the state \rp, 77) does not represent nature. It may then be considered as a 
superposition of possible states of nature, but there is no guideline how to determine the states in the superposition. 

Moreover, in the usual approach, (tp, t]d\^(x.)^/(x.')\^P, T)r>) {vd represent the conformal time of decoupling) can be 
readily shown to be rotationally invariant, a function of x • x'. Therefore, again, \ip, 77) does not represent the state of 
our universe, but at best is some superposition of possible states. 

So, as we have emphasized, the homogeneous and isotropic initial state and dynamics does not explain the observed 
inhomogeneity and anisotropy. 

Therefore, as previously mentioned, following [Tl|[l| we take a different approach, utilizing the semi-classical descrip- 



tion of gravitation 34] 35] , where gravity is treated classically while other fields are treated in the standard quantum 



field theory (in curved space-time) fashion. The classical gravity and the quantum fields are thus related by 

G ab =8irG(f ab ). (42) 

There is an immediate objection to semi-classical gravity. Suppose a quantum experiment is performed with two 
possible macroscopic outcomes, a large object being put in one or another place. Using the Schrodingcr equation 
to describe this, including the apparatus, the resulting state vector describes a superposition of these two outcomes. 
Then, (Tab) is large in two places and so, according to semi-classical gravity, the gravitational field acts as if there 
were sources in two places. Such an experiment was actually performed 36] and, as expected, the gravitational field 



The rcd/bluc shift contribution is ^~ = -^yM -. a >1> where 9go'°' represent the " time-time" metric components at the events of 
emission and observation respectively. For details about the -2/3 contribution, see l23j p. 139, or |24H . 
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(as measured by a Cavendish balance) saw only a source in one place. However, this objection no longer obtains if the 
Schrodingcr equation is modified, a la CSL, to include collapse, since then the state vector rapidly ends up describing 
the object in one place onhT^l. 

However, the resolution of this problem by invoking collapse brings on another problem. As is well known, intro- 
ducing collapse violates the conservation of the energy-momentum tensor. If the divergence of the energy-momentum 
tensor does not vanish, and it is equated to the Einstein tensor, then of course the latter's divergence does not vanish 
either. 

Here we shall take the view that Einstein's equations are an emergent, approximate description of the collective 
behavior of the fundamental degrees of freedom of quantum space-time, and as such those equations will not hold 
under all circumstances. It has been argued in [l3j that this should be considered in analogy with the breakdown 
of the Navier-Stokes characterization of a fluid. This can be expected to occur, not only for phenomena at scales 
smaller than the mean intermolecular distance of the fluid constituents, but also when there are important energy 
fluxes between the micro and macroscopic degrees of freedom, such as when a part of the fluid undergoes a phase 
transition. 

In the same manner, the collapse should be thought of as accompanied by a back reaction in the fundamental 
quantum gravity degrees of freedom, which are not fully represented in the metric characterization. Then, a more 
precise description would include a compensating term appearing in the Einstein's equation. Another approach, 
discussed for instance in (38j , involves assigning energy-momentum to the stochastic field driving the CSL dynamics 
such that the stress tensor then does have vanishing divergence. We shall not explore this issue any further in the 
present work. 

It follows from Eqs.(gT])and Eq.([I2]) that 

- kH( Vl k) = 47rG$,(»jwy(M)) = A * G< &w (43) 

((7r(k, 77)) = (i(j,r]\Tr(k)\t/j,r])). When inflation starts at time 77 = — T, it is supposed that the state is described by 
the Bunch-Davies vacuum, so (tp, — T|7r(k)|^, — T) = 0, and the space-time is homogeneous and isotropic. While this 
would remain the case were there only hamiltonian dynamics, the addition of CSL dynamics causes (-7r(k, 77)) 7^ 
thereafter. The CSL hamiltonian depends upon a classical random function of time w(t) (more precisely, one w(t) 
for each momentum mode). Each set of such 7«(t)'s gives rise to a different possible inhomogeneous and anisotropic 
universe. 



V. THE OBSERVATIONAL QUANTITIES 

The quantity that is measured is AT(9,<p)/T, which is a function of the coordinates on the celestial two-sphere. 
This data is expressed in terms of spherical harmonics as 

= Y,- lm Y lm {6M, m m = J d^^My^eM- (44) 

Ira 

**We emphasize, as already discussed, that we are factoring out the late time physics, so we discuss 'observations' 
were the transfer functions equal to constants. 

The experimental results are usually expressed in terms of the quantity 



m 

What is 'observed,' then, is that the quantity 



Cl z 21 



OB l = 1(1 + 1)(2Z + l)- 1 £ K„| 2 = 1(1 + l)d (46) 

in 

is essentially independent of I. This 'scale invariant' (the reason for the name and the arcane dependence on I shall 
be given subsequently), or 'Harrison-Zcl'dovich' spectrum, is what must be accounted for by the theory 



For more discussion about the applicability of semiclassical gravity to the problem at hand see [l3ll where the first steps of a formalism 
capable of incorporating collapse of the wave function at the semiclassical level was developed. 
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Our approach produces explicit expressions for the quantities that are most directly extracted from the data. Using 
Eqs. ([gg j) .(|ig ]) , we obtain 



-c / d d /ce Ikx — (7r(k,7j)), where c = ruw£_ (47) 

7 fc 3a 



r 



Here, x is the coordinate of the point on the intersection of our past light cone with what will eventually become the 
last scattering surface in the direction on the sky specified by 8, (p. Then, according to Eq. (|44"|) . 



a lm = cj #toY? m (6,<p) J d 3 fee ik ' x -i<7r(M)>. (48) 

Thus, aim depends upon (7r(k, 77)), a well-defined stochastic quantity in our treatment. (This differs from the standard 
treatment, where no comparable expression can be given.) The stochasticity occurs because the collapse theory gives 
an ensemble of possible universes (only one of which is actually realized) and the associated probabilities of realization. 
We shall see that there is not a large deviation from the mean, so we may consider that our universe is typical. 
Continuing, it follows from Eq. (|4"5)l and the well-known expansion e lkx = 47r^ ; m i l ji(kr)Yi m (8,(p)Y l l n (k) that 



1 



a lm = i l 4nc J ^kjiikR^Y^k)^^,^). (49) 

Here, Rd is the co-moving radius of the last scattering sphere, so x = Rd(shi(0) sin(i^), sin(0) cos(</>), cos(0)), and k 
is the unit vector in the direction k = kk. Therefore, 

\a lm \ 2 = (Aire) 2 J d 3 M 3 fc'iK^)iKfc'i2 D )^ r „(fc)y4(^0^((*(k,r?)>(^(k',r ? ))*). (50) 

We may consider that the average over the ensemble of possible universes fairly reflects the value obtained in our own 
universe. As shall be seen, the form of our expression for the ensemble average at the end of the inflationary period 
is ((7r(k,77))(7r(k',77))*) = f(k)S(k - k'). Then, 

r 1 r°° 1 

\ai m \ 2 = (4^c) 2 J d 3 kj 2 (kR D )\Y lm (k)\ 2 -^f(k) = (4ttc) 2 J dk 3l {kR D f -f(k). (51) 

Now, we note that if f(k) = ak where a is a constant, the result becomes independent of Rd- That is what is 
referred to as a scale invariant spectrum. In that case, 

r°° 1 1 

\a lm \ 2 = (^c) 2 aj^ dx J?( x )- = (^ c ) 2a 2l{l + 1 y ( 52 ) 
Thus, our estimate for the quantity that is usually the focus of the analysis is 



Cf = ^Ekf = N^. (53) 

Thus we obtain 1(1 + l)Ci is a constant, independent of I, in agreement with 'observation,' as mentioned following 
Eq (|46[) . However, as we have seen, that only occurs if 



(ir(k)} 2 - k. (54) 

So, we now turn to add CSL dynamics to the dynamics discussed in Section IIIIDI (which is governed by the 
hamiltonian H/., Eq.(|28p). to see whether, under the combined dynamics, the ensemble of possible universes can 
satisfy (Eq.([5gj)). That requires (P) 2 - k. 
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VI. CSL 



We shall be using just the simplest form of CSL, which describes collapse toward one or another eigenstates of an 
operator A with rate ~ A. 

As we have seen, the relevant operators on which wc should focus our attention are the 7r(k, rf) ~ P, and their 
expectation values (7r(k, 77)) ~ (P) 7^ 0. Wc may call P our 'focus' operator to differentiate it from the 'collapse 
generating' operator A. (Note that we could choose A = P, or make another choice for A.) 

There are two equations we must consider. 

The first is a modified Schrodingcr equation, whose solution is: 

IV,*) = re-^^+Awo-^ip] |^ 0) _ (55) 

(7~ is the time-ordering operator). w(t) is a random classical function of time, of white noise type, whose probability 
is given by the second equation, the Probability Rule: 

PDw(t)^^M,t)f[^^=. (56) 

The state vector norm evolves dynamically (does not equal 1), so Eq. (|5"6f says that the state vectors with largest 
norm are most probable. That the the total probability is 1 can be seen from 

/ PDw{t) = [ Dw{t-dt) H dw{t) (V>, t - dt\e~ dt ' e - fj dt' [iH + ^Mt')-2XA]'] ^ f _ dt) 

J J J-oo y'2'KX/dt 

= j Dw {t-dt) J^^^= {ij ,t-dt\e~^'^ 2 \^t~dt) 

= J Dw{t-dt){tp,t-dt\tp,t-dt) = ... = (V>,0|V>,0) = 1. (57) 

To see how the dynamics collapses to eigenstates \a n ) of A (assuming H = 0), write \ip, 0) = Y^n=i c n\ a n) so, 
according to Eqs. (j55|) . f[56 ]l . 

TV N 

|V,*) =e~^fo w2 V)J2cnK)e B ® a "- xta2 % P = e-A^- 2 (0^| c 2| e 2S( t )a„-2Ata^ (5g) 

n—1 n—1 

where B(t) = J * dt'w{t'). Writing w(U) = B(ti + dt) — B(U), so Y[dw(U) = Y[dB(ti) : we can integrate P over all 
B(ti) except B(t), obtaining the result 

P'(B(t))dB(t) = V | c „| 2 ^£e"^[ s W- 2Ata "] 2 . (59) 

n—1 v 

According to Eq. (j59")) , the probability is the sum of gaussians, each drifting by an amount ~ a n t, but of width 
~ \f~Xt. Therefore, after a while, they evolve into essentially separate gaussians. Then, there are ranges of B(t) which 
correspond to each possible outcome. If —Ky/Xt < Bit) — 2\ta„ < Ky/Xt, (K > 1 is some suitably large number), 
the associated probability integrated over this range of B(t) is essentially |c„| 2 , and the state vector given by Eq. (|55|) 
is essentially \ip,t) ~ \a n ). 

It should be emphasized that, when H ^ 0, the hamiltonian dynamics interferes with the collapse dynamics. In some 
cases, collapse nonetheless takes place. **In some cases, a kind of stasis or equilibrium between the two competing 
dynamics is reached. In other cases, the unitary and non-unitary dynamics interfere with each other in interesting 
ways. 

It is useful to have an expression for the density matrix which describes the ensemble of evolutions. This is obtained 
from Eq.(l55l): 

' Dw(t)Te- /o* dt ' [i£+&MO-2A4»] ^ Q) ^ Q|e - ft it ' [-iH + ^ [w(t ')-2XAf] 
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where the subscripts L and R mean that the associated operators are to be put to the left or right of p(0), and the 
T reverse-time-orders operators to the right of p(0). The evolution equation for the density matrix is therefore 

±p(t)=-i[H,p(t)}-^[A,[A,p(t)]]- (61) 
It follows that the ensemble expectation value of an opcrator(O) = TrOp(i) satisfies 

|<3) = -i\6Jf\ ±[A,[A,d}}. (62) 



VII. APPLICATION OF CSL 

CSL dynamics tries to collapse state vectors toward eigenstates of A. It gives an ensemble of different evolutions 
of the state vector, each characterized by a different w(t). It is to be applied to the modes described by the focus 
operator P, the operator X, and hamiltonian given in Eq. (|2"5]) . According to Eq. tfM]) . using the relations which 
gave us Eq. (j3"7) . 



(7T(k))(7r(k')>* = ((7T S (k) +^(k)))((7T S (k') - Z7T A (k'))> = (fa(k)} (tT S (k')) + (* A (k)> <7TA (*)) 



= 2(-=) 2 ( 5 kk , = (P) 2 <5(k - k'), (63) 



what we need to find is the ensemble average (P) 2 , and determine under what circumstances, if any, this is ~ k. We 
must therefore choose a collapse generating operator A. We shall consider the simplest possibilities in this paper, 
A = X and A = P. These correspond to the basic inflaton perturbation operators 8ip(x), Stp'fa). 

Because, in these cases, Eq. (|62[) is quadratic in the operators of the hamiltonian and of the collapse hamiltonian, 
one can readily obtain a set of coupled equations for the ensemble average of any power of operators, just as was done 
in Section Imll] But, 



<<M|P|V,7?> 2 _ f n i A^,v\P\i>,v) 2 



<P> 2 = / P^) 7.". y ' = / Dw(ri) v ' |y ' " (64) 

can not be obtained in this way, and is not so easy to calculate directly. However, for this problem, there is a 
relationship between (P) 2 and (P 2 ) which allows us to obtain the former more easily. 

Because the initial state is a gaussian and the hamiltonian and collapse hamiltonian are quadratic in X, P, the 
form of the state vector in the momentum basis at any time is 

<p|V, v) = e - AMp2+BMp+CM . (65) 

The initial conditions are A(— T~) = l/2fc, B(—T) = C(— T) = 0. (By solving the Schrodingcr equation with this 
ansatz, which we shall eventually do, one finds that A is independent of w(t), B is linear in w(t) and C is quadratic 
in w(t).) Therefore, the momentum matrix element is 



1 (B + B*) (B+B'f , r , r ^ 



2y / (A + A*) (A + A*) 

We also note that the state vector norm is 

ty^T,) = f dpe -(A+A-) P 2 +(B+B*)p+(C+C*) 

Therefore, by Eq.(|64| 



( S -\- S * ) ^ 
j 4(A + A*) p(C+C ) 



V(a + jF) 



(p) 



Dw(rj) 



1 



(B + B*) 
2(A + A*) 



gi(A + A*) g(C + C*) ^ 



(66) 



(67) 



(68) 
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Now, (P 2 ) can be written as: 



(P 2 



Dw I dpp 2 e^ A+A '^^ B+B '^ c+c ' ) 



Dw- 



1 



2(A + A*) 12{A + A*) 



(B + B*) 



g 4(A + A') giC + C) 



2(A + A*) 



(Py 



(69) 



where the last step follows from (A + A*) being independent of w(t), from the integral of Ea.(|67p being 1 as it is the 
probability of all possible u>(i)'s (Eq.(j57j)), and from Eq.([68]). 
To summarize, 



(P) 2 = (P 2 ) 



1 



2(A + A*) 



(70) 



i.e., [2(A + A*)} 1 is the standard deviation of the squared momentum. It is also the width of every packet in 
momentum space. 

Thus, to calculate (P) 2 , we shall find the second term on the right hand side of Eg . ([70)1 from the Schrodinger 
equation and we shall find the first term by using the density matrix . 

VIII. P AS GENERATOR OF COLLAPSE 
A. Use of Schrodinger equation. 

The Schrodinger equation is the time derivative of Eq. (|55p . In the momentum representation, with A = P, it is 



dr] 



p 2 - -lp— + —p) - k 2 — 
rj dp dp dp 2 



<p|^M> 



w 2 (v) ~ w(r])p + \p 2 }{p\il),r]) , 



We note that A is a dimensionless number. The equation A satisfies is 

^A=[- + \]--A~2ik 2 A 2 . 
dr) 2 T) 

This Ricatti equation is solved by writing A = Z/[2ik 2 Z]. Putting this into Eq. ([T2)) we get: 

i]Z = 2ik 2 [^ + \]r)Z -2Z. 



Defining a = fcyl — 2i\, the two solutions are 



Therefore, 



A(v) 



2ik 2 rj 



— cos ar\ and — sin an. 



cos ar\ — arj sin arj + C (arj cos ar\ — sin arj) 



cos arj + C sin arj 



(71) 



(72) 



(73) 



(74) 



(75) 



The constant C is determined by the initial condition A(—T) = 1 /2k: 



C 



(1 — ikT) cosaT+ aTsinaT 
(1 — ikT) sin al~ — aTcos ofl~ 



(76) 



16 



Putting C into Eq.flTSJl yields 



2k 2 Tj 2ik 2 



(1 — ikT) cosa(?7 + T) + aTsma(r] + T) 
(1 — ikT) sin a{rj + T) — oT cos a{q + T) 



2k 2 ?] 



a 

2k 2 ' 



(77) 



The last step has utilized the approximations fcT >> 1, << T, a ~ k — iXk for A << 1, yet AfcT » 1 (to be 
justified in section IX D[) so cosa(?y + T) ~ isma(rj + T). Therefore, we have the result, at rj = — r, that ** 



2(A + A*) a + a* y/\ - 2i\ + y/l + 2i\ yfiy/i + y/T+W 



(78) 



B. Calculation of (P 2 ) 

For this case, Eq. ([62l) becomes 



-(0) = -i[0,H]--[P,[P,0]. (79) 

Referring to Scctior llllDl where we considered the hamiltonian dynamics alone, the first order Eqs. (|31[) arc unchanged, 
so here too (X) = (P) = 0. 

The second order equations Eqs. (pt5|) for Q = (X 2 ), R = (P 2 ), S = (XP + PX) are unchanged except for the first: 

Q = S- — + X,R= -k 2 S + — , S = 2[R - k 2 Q}. (80) 
7] rj 

The general solution is therefore the sum of the three solutions Eqs. ([M|) to the homogeneous equations added to an 
inhomogeneous solution: 

Q = Xrj/2, R = Xi]k 2 /2, S = A/2. (81) 
For example, the equation which replaces Eq. (|34b[) is 

\k 2 n 

R = Cie 2ikr, + C2e -2ik V +Ca + ^_H_ ( 82 ) 

The constants are determined by the conditions at t = — T, which are Q = l/2fc, R = k/2, S = 0, . Assuming 
fcT >> 1, it follows from the modified Eqs. tjM)) : 

— = C\ + C*2 + C3 — 

1 _ r 1 r 1 r 1 AT 

{)=-2iC 1 \ + 2iC 2 \ + ^-. (83) 



The solution is 



k k 2 

Putting Ea.([55p into Eq.([52]l, setting rj = -r, and using fcr << 1, fcT >> 1 we get: 

Afc 2 T fc , , 

R = — +2- ^ 

When AfcT >> 1 it is clear that the first is the dominant term, but we include the second term to enable consideration 
of the case A«0. 
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Therefore, by Eq.flTO]), using the results (JTSJ) and ([85]). we obtain 

/p\2 _ ^k 2 T _j_ & & 

2 2 V^v 7 ! + VTTiA 2 



(86) 



If we set A = (turn off CSL), we have the standard quantum mechanics result (P) 2 = 0. 

We see that agreement with the observed scale- invariant spectrum, (P) 2 ~ k, can be achieved if we assume the first 
term is dominant and we also set 



A = A/fc. 

We note that this replaces the dimcnsionless collapse rate parameter A with parameter A of dimension time -1 . 
In that case we obtain: 



(87) 



(P) 



XkT k 



k 



2 2 



V^/l + 4(A/fc) 2 



It is also worth noting that, if A << 1 but still AT >> 1, then 

{[P-{PW i 



« 1, 



(P 2 ) AT 

and the universes in the ensemble do not deviate much from each other. 



IX. X AS GENERATOR OF COLLAPSE 



This proceeds in a manner parallel to the previous section. 



A. Use of Schrodinger equation. 

In the position representation, with A = X, the Schrodinger equation is 
d 



We note that A has dimensions time 2 . The wave function in the position representation is {x\ip,ri) = exp[— A'x 2 - 
B'x + C], where A' — 1/4A (A is the coefficient of p 2 in the exponent of the Fourier transform of (x\ip, ^satisfies 



±A' = [— + A] + —A' - 2iA' 2 . 
w) 2 r\ 

This Ricatti equation is solved by writing A' = Z/[2iZ]. Putting this into Eq.([W|) we get: 

r)Z = -/3 2 i]Z + 2Z. 



Defining (3 = \/k 2 — 2i\, the two solutions are 



Using Z = /3 2 ryexpi/3ry, 



A! {n) 



2i 



e 2iPv + c 



s 2i ^(l- ipn) + C{l + i/3ri) 



(90) 



(91) 



(92) 



(93) 
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The constant C is determined by the initial condition A'{— T) = k/2: 



C 



D -2if3T 



/3 2 T - k/3T + ik 
/3 2 T + k/3T + ik ■ 



Putting C into Eq.® yields 



A! {n) 



2i 



e 2^(T+r,) [ / g2 T + k/3T + ^ _ [ / g2 7 - _ kj3T + - fc ] 

e 2ip(T + V ) [p2j- + k pq- + ik] (i _ ipr,) - [p 2 T - k/3T + ik](l + i/3rj] 



'V 



2(i + Prj) ' 



where the last step has utilized fcT >> 1, and exp 2(— Im/3)T >> 1. 
Using A = 1/4,4', we find 



1 



W = . (l + 4(A/fc 2 ) 2 ) 

2(A + A*) Rc(A') [ 1 > F(\/k 2 ) + 2(\/k 2 ) 2 F- 1 (\/k 2 ) - 2{\/k 2 )(kr ] )- 1 ' 



(94) 



(95) 



(96) 



where F(x) = -^Vl + VTTix 2 (T(0) = 1) 



B. Calculation of (P 2 ) 

For this case, Eq. (f62"l) becomes 



-(G) = -i[G,H]-'-[X,[X,G]. (97) 

Referring to Scctior llllDl where we considered the hamiltonian dynamics alone, the first order Eqs. (|31[) arc unchanged, 
so here too (X) = (P) = 0. 

The second order equations Eqs. ([33|) for Q = (X 2 ), R = (P 2 ), S = (XP + PX) are unchanged except for the 
second: 

Q = S - — , R = -k 2 S + — + A, S = 2[R - k 2 Q}. (98) 
77 77 



The general solution is therefore the sum of the three solutions Eqs. (|34|) to the homogeneous equations added to an 
inhomogeneous solution: 

Q = Xr//2k 2 , R = A77/2, S = 3A/2fc 2 . (99) 
For example, the equation which replaces Eq. (|34b[) is 

R = Cie 2ifc " + C 2 e~ 2lk7 > + C 3 + ^. (100) 

The constants are determined by the conditions at t = — T, which are Q = l/2fc, R = k/2, 5 = 0,. Assuming 
fcT >> 1, it follows from the modified Eqs. (|3"4")) : 

= -2^1 + 2^1 +3 A (101) 



The solution is 



« = f + f. (io2) 



Putting Eq. (|102[) into Eq. (|100[) , setting r\ = — r, and using fcr << 1, AT >> k we get: 



«=f4 (103) 
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Therefore, by Eq. ([70)) . using the results (|95|) and (|103[) . we obtain 



AT Vfc 4 + 4A 2 XT k (l+4(X/k 2 ) 2 ) 



{P)2 2 V^ 2 - 2iX + sjk 2 + 2iX 2 + 2 (fc/2) F(A/fc 2 ) + 2(A/P) 2 F- 1 (A/fc 2 )-2(AA 2 )(fc77)- 1 - (1 ° 4) 



Once more, if we turn off CSL, we find (P) 2 = 0. 

We see that agreement with the observed scale-invariant spectrum, (P) 2 ~ k, can be achieved if we assume that 
the first term dominates, and if we set 

A = Xk. (105) 

We note that this replaces the collapse rate parameter A of dimension time -2 with the parameter A of dimension 
time -1 . In that case we obtain: 

AfcT k (1 + 4(A/A:) 2 



{P) 2 = ^ + ^\ l = ~ ]• (106) 

2 2 L F(A/fc)+2(A/fc) 2 F -1 (A/fc)-2(A/fc)(fc? ? )- 1 



The validity of our approximations shall be discussed in Section [XDI 



X. PHYSICAL QUANTITIES 



As we have shown, our theory can describe the Harrison-Zcl'dovich scale-invariant spectrum. Moreover, it can be 
used to find expressions for various physical quantities, such as a; m , AT(n), ^(x) (and their probabilities of taking 
on various values), which is not possible with the usual approach. This we shall now demonstrate, using collapse 
generator A = P discussed in Section [Villi 

First, we need the expression for (?/>, r)\P\^), rf) / (if), T)\i/>, rf) = (P). To get this, we return to the Schrodinger equation 
Eq. dTTI) for (p\ip, rf) = exp[— Ap 2 + Bp + C] and, with use of Eq.([T7j). obtain the equation for B: 



—B = --B - 2ik 2 AB 
dr\ rj 



—B- 2ik 2 
V 



i a 
2k 2 r] + 2k 2 



Thus, according to Eas.(|B5 p .l|57 )l . and 1(1070 . 



(P)(V) = 



B + B* 
2{A + A*) 



R 



= — I drfw{rf)e 



B + w = -iaB + w or B{t) = / dr/w(r/)e - * a( ' 7- ' 7 ' ) . 

(107) 

(108) 



'\ p -S(ri-ri') 



cos R{q — t]') 



T 



where we have written 



= fcVl - 2iX = R-iS, so R 2 -S 2 = k 2 , RS = Xk 2 , so R = fcP(A), S = 



Xk 
W)' 



(109) 



Now we can express physical quantities in terms of (P). As was done in obtaining (|63[) . and using (|108[) we write 



(7r(k,r/)) = (71-5 (k, t?)) +i(n A (Kil)} 



2Vdk 



2R 



drf[w s (k, rf) + iw A Qt,V , )]e~ S(ri ~ ri ' ) cosR{?] - rf), 



-T 



(110) 

where we have introduced white noise functions w(k, rf). These are only defined in the upper half k-planc. However, 
the fourier transform of 7r(k) is 7r(x) which is real. This implies wg(— k) = ws(k) and wa(— k) = — -^(k), so Eq. (|110[) 
holds for all k. 

To go along with the expression (|108[) for (P), we must have the probability of w(rj). To get the probability, we 
could go back to Schrodingcr's equation Eq.([7T|) and calculate C. However, since the probability is (ip, t)\i/j, rf), it is 
easier to use Eq.([7T|) to get 

r?|V, rf) = ±^{i>, ryl-0, rf) + Mv)(1>, l\P\1>, v) ~ 2A(^, ry|P 2 |^, V ) 



+ 2w(r 1 )(P) 2A(P 2 )] fanW,*) = { ±[w{ V ) 2A(P)] 2 - ^j^r^r?) ( 



111) 
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where the last step follows since Eq |70|) holds without the ensemble average. Thus, we obtain the probability density 

P{w) = (iP,t\xIj,t) = e -^^T d nMv)-2HP)(n)]\ (112 ) 

where (P){rj) is given by Eq. (|108p (the factor exp — \k{r+T)/ F{\) has been absorbed in the normalization appropriate 
to Dw). 

To go to the continuum case, replace J dr\ by J drjdk and w(rf) by u>(k, rj) in Eq. (|112| ). 

In order to do calculations with the probability (|112p , it is convenient to define a new random variable: 

rv 

v(r)) = w(r)) - 2A(P)(r?) = w(r)) - 2S J dr/w^e- 8 **-" ] cosR{r] - rf). (113) 

The Jacobian determinant of the transformation from variables w(rj) to 17(77) is 1 (essentially, the matrix transformation 
has all diagonal elements =1, and zeros to the right of the diagonal). The probability density of v{rf) is very simple, 

P(v) = e-AZ-T^Mt?)^ ^) = 0, v{r,)v{r,') = \5(r]-r)'). (114) 

However, since the physical quantities are expressed in terms of w, we need to invert Eq. (|113j ) to obtain the expression 
for w in terms of v. This is done in Appendix \A\ with the result: 

w{rj) = v(?]) H — — / drj'[Ssmk(7] — rf) + kcosk{j] — r]')]v {rf). (115) 
k J-t 

We shall provide a few examples of the use of this formalism. We shall show, using it, that we obtain the result (|55|) 

for (-P) 2 . We shall calculate the probability distribution of the temperature fluctuations (which shall prove to be a 
gaussian) as well as the correlation function of the temperature fluctuations. Finally, we exhibit the expression for 
ai m , and the ensemble average of |a; m | 2 . 



A. Calculation of (P) 2 — Again. 



As a consistency check, we employ a different way of calculating (P) 2 than was done in Section lVIIII Using Eqs. (|108[) 
and (|115[) , and taking r\ = 0, we write 



k 2 r° 

(P) = — drfe S71 cos R{ri') 
R J-T 



2S f 71 

v iv') H — r~ / drji[S sin k{r]' ~ rji) + k cos k{r]' — rji)]v{rji) 



(116) 



We immediately note that (P) = 0, since v{rf) = 0. 

The order of integration of the double integral can next be exchanged, and the integral over 77' performed. There 
is a term which cancels the single integral in Eq. (|116[ ). and the result is 



(P) = — J drjiv{rji) [k cos krji — S sin krjxj . 



-T 



Therefore, with use of Eq. f|114[) , and neglecting terms small compared to kT, we obtain 



(P) 2 = — J d m [kcosk Vl - Ssmk m ] « [S 2 + k 2 ] = — ^— . 



(117) 



(118) 



This is the same result as in Eq. (|86 



B. Temperature Fluctuation 

The temperature fluctuation at the end of inflation is given by Eq. (|4T|) with 77 = — r. (Since kr << 1, we shall 
replace r by 0.) With use of Eq. (|110p ). we therefore have 

AT f d\c ( die k 2 P° , 

^jr= c -^e 4MDk "(7r(k,0)) = c j _ e ^ kfi — / drf[w s {k, rf) + iw A (k, r)')]e Sr > cosfiry' 



r 



r dk k 2 f" 

/ -z— I drfe Sr] cosR{r]')[cos{kR D k ■ h)ws{k, rf) — sm{kRok ■ n)wA{k, rf)} (H9) 
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To conveniently calculate probabilities, we need to replace the w's by v's, using Eq. (|115[) . Just as in the previous 
section, we may then exchange the order of the double integral, obtaining 



AT(n) 



f dk k f° 

I ~~i2~5 I ^i[cos(fei?£)k • n)us(k, — sin(fci?£)k • n)u^(k, r?i)] [fccos(fc?7i) — /S'sin(fcryi)] . (120) 
J + k R J -T 



1. Probability Distribution of the Temperature Fluctuations. 
The probability that AT(n)/T = C is given by 

P{G)= [ Dve- f r d vU dk ^l v sM+ v2 ^S(AT(n)/T-C). (121) 



Writing Eq.([T2DJ) as 

AT(n) 



dm f dk[f s {k 7m )vs{k, m ) + f A (k, m )v A (k 7m )}, (122) 
-T J + 



T 

Eq. (jl21|) becomes 

P(C) = j Dve^^^f+^^^^+^^—J 00 rf a;e ^(/-r d7 'i/+ dk [/s(k,»?i)^s(k,7 ;i )+/ A (k > 7 ;i )«A(k^i)]-C) 

= — f due~ iu)C e~ u2 $- r dm A- dki T 1 [-' ; I( k ' ? ' 1 )+-' ? -4( k ^i)] _ e _ ^^T d "i /+ rfkA(fcu/g(k,, J1 ) + /;t(k,„ 1 )] 
Thus, we see that the probability is gaussian, with variance 



AT(n) 



dk\(k)j—^J ^d m [kcos(kr h ) - Ssm(km)] « c 2 j dk\(k)^ = tic 2 XT ^ — . (124) 



r 

We have set A(fc) = A/A; in the observed range of k, 10~ 3 Mpc~ 1 < k < W 2 Mpc~ 1 , in order to obtain agreement 
with the Harrison-Zel'dovich spectrum. However, [AT(n)/T] 2 is a measured quantity so it seems reasonable to expect 
that the integral in Eq. (|124j) should converge, and it doesn't. 

One might think the issue of the divergence of the above integral could be resolved simply by taking into account 
the fact that the temperature fluctuations seen in the CMB, are the result of modifications by "late time"El] processes 
of the primordial seeds of cosmic structure we have calculated. The fact however, is that such modification^] occur 
after the reheating stage, and are therefore of no help at the time just prior to reheating, which is the era to which 
our calculations refecJ- 

If there were no modification in this estimate for the inflationary regime itself, the resulting ensuing fluctuations 
would be predicted to be of arbitrarily large magnitude and a perturbative treatment would, of course, be invalid. 

The issue is intimately tied with the universal scale invariance of the H-Z spectrum and is thus a problem that, 
although often unrecognized, afflicts the standard treatment of inflationary perturbations. Let us consider for instance 
the treatment presented in [24|. 

Consider the expression for the Newtonian potential at a point x. To do this one uses Eq.( 8.5) of that book and 
writes: 



11 That refers to the regime well after inflation has ended. 

12 These are due to well understood physics connected with the behavior of the plasmas of ordinary matter (quark -gluon plasma before 
hadronization forms the nucleons, nucleon-electron- photon plasma before nucleosynthesis, and Hydrogen-Helium plasma after nucle- 
osynthesis) codified in the so-called transfer functions which describe how each mechanism alters the temperature structure. For example, 
transfer functions are known to exhibit a damping effect known as "Silk Damping" at short distances (large k) due to viscosity/photon 
diffusion. 

13 Strictly speaking at the end of inflation the conditions are still far from thermal equilibrium, there is therefore no temperature and 
in fact no radiation. It is only after the reheating that the decay of the inflation is supposed to lead the system towards a state of 
"thermal equilibrium" . However the quantities we are computing have their counterpart as the inflaton-energy-density fluctuations and 
Newtonian potential fluctuations. 
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The ensemble average of the square is then, 



$ 2 (x) 



(27T) 3 



1 



/ 



dkdk'$ k $*,e t ( k - k ') x . 



(126) 



Next one uses the definition 8.11 of [24| which indicates that we might write the ensemble average as, 



$ k $£, =6(k-k / )2Tr 2 5l(k)/k 3 



(127) 



where in their notation 5%{k) is the power spectrum (for the gravitational potential $) and is taken to be a function 
only of k = ||k||. Thus, substituting in 11261 we find: 



Now we take eq 8.103 of [24[ which indicates that the spectrum is scale independent. This requires 5\(k) to be 
constant, independent of k and, in fact, given by 4(e + p)/C s where e and p are the energy density and pressure of 
the background inflaton field respectively, and where C s is the effective "speed of sound" (see cq. 8.50 of [24[) which 
for a canonical scalar field is just the sped of light. We thus obtain: 



which diverges just as our estimate of the temperature fluctuation in 11241 

So, it seems that the problem we face here is not intrinsic to the CSL theory or, more generally, to the hypothesis 
that some kind of collapse of the wave function plays an important role in inflationary cosmology. What the analysis 
in the present paper does is show the problem more explicitly. After all what the collapse does, roughly speaking, is 
provide a mechanism to convert quantum mechanical uncertainties already there in the initial state into a range of 
actual values at a later time. 

Of course, if the problem lies beyond the collapse approach and seems intrinsic to the inflationary proposal for 
generation of primordial cosmological inhomogeneities and anisotropies, that does not mean it can be ignored. One 
approach might be to consider a similar problem that appears in any treatment of quantum fields: the uncertainty 
in the value of a field in its vacuum state at any point in Minkowski space-time is infinite. This is often taken to 
indicate that one cannot consider such a quantity as the field at a point and must, rather, focus attention on things 
like smeared fields over space time regions. After all, quantum fields are distribution-valued operators rather than 
ordinary operators on Hilbcrt space. In the inflationary context however it is rather unclear what would be the 
appropriate smearing one should consider. 

For our problem, we might rather consider the possibility that (|124p be made finite as the result of a more compli- 
cated behavior of X(k) ^ A/fc as k — > and k — > oo. For example, we might set 



for ki >> k\ and the interval (fe, k\) very much wider than the range of k appropriate to the CMB. That would in 
effect make a prediction, that the H-Z scale invariant spectrum not hold outside a particular range of k. 

On the other hand this issue might be resolved by things unrelated to the CSL collapse rate and therefore using 
the above argument to modify it might be premature. One possibility that would certainly help in this regard is the 
known fact that, during the slow roll inflation, the expansion rate is not exactly de-Sitter like. This leads, in the 
usual analysis, to a modification of the H-Z spectrum corresponding to a factor where the "spectral index" 

n s is known both theoretically and empirically to be smaller than l[37}. This would remove the divergence coming 
from the k — > oo behavior. Regarding the divergence arising from the k — > behavior it seems that the most natural 
resolution would come from noting that the region that undergoes inflation would in all likelihood have started as a 
a very small patch, which despite the nearly exponential expansion resulting from inflation will remain finite. Thus 
the range of physically relevant values of k would be bounded from below by some kfi n i te > 0. 




(128) 




(129) 




(130) 
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2. Correlation of the Temperature Fluctuation. 
The correlation function of the temperature fluctuation can be found from Eq. (|120[) 

A ^ (n) AT ^ =c 2 J dkX(k)j^cos[kR D k-(n-n')} J° d Vl [k co S (k m ) - S sm(k m )] 2 

pOO 

= 4(nc) 2 T / dk\(k)Y'j?(kR D )Y lm (n)Y l * m (ti) 

70 lm 
poo 

= nc 2 r dk\{k)Y j (2l + l)jf(kR D )Pi(n-n / ). (131) 
Jo , 

(Note, XM 2 ' + l)i? fa) = !> so Eg.(fl3T]) agrees with Ea.([T2ijl when h = h' .) 

C. ai m 

From the definition of a; rn and Eq. (|120[) we can exhibit the expression for ai m : 



ah 



J dk-^ j ^ dm ^(2Z + l)ji{kR D )Pi{k ■ n')[S(l)v 3 (k, m ) + iO(l)v A (k, m )} [kcos{k m ) -Ssm(k m )] , (132) 



where 8(1) is 1 if I is even and if I is odd, and 0(1) is if I is even and 1 if Z is odd. 
However, it is easiest to use the second of Eqs. (|131[) to find 



|a im | 2 = (2ncj 2 T j dk\(k)jf(kR D ) = (2^) 2 AT J T jf(kR D ) = ■ 



" dfc -2/7 d \ (27rc) 2 AT 



which is to be compared with Eqs. (f5~T|) . (|52"j) . In Eq. (|133p we have set A(fc) = A/fc, although we earlier specified that 
it might behave differently as k — > and k — > oo. This different behavior has been ignored in calculating the integral 
in (|133[) . which converges perfectly well without it. Thus, such possibly different behavior is now constrained to have 
a negligible effect on the integral in Eq. (|133|) . 

This concludes our demonstration that the dynamics produces agreement with the observed Harrison-Zcl'dovich 
spectrum. 

D. Estimates 

It is necessary to make some order of magnitude estimates. We start from the fact that the temperature fluctuations 
in the CMB are ^f- = 1/3^ ~ 10~ 5 . We can now use that together with our other results to check the self-consistency 
of assumptions made in earlier sections. We note that we took the dominant term in Eq. ([88|) to be 

XkT 



(P) 2 * -y-- (134) 
This led to the result (TT!M|) . 

2 ~\TI (135) 



r AT" 






I T . 


~ 4 


3a . 



where X = 1/A J dk\(k) is at least as large as J^ °- 3 dk/k w 11.5. 

Using the definition % = a'{rj)/a(ji) in the second of Eqs.([5]) and taking the derivative with respect to 77 we find 

mm! = 47rG(2^Vo + 4aa'V[fo] + l^d^V^oWv). (136) 

Substituting here the expression for from the first of Eqs.© gives 

mW = 4irG(-4H{ct) ) 2 + Aaa'V[4> ]) = l^GU(-(<j)'uj 2 + a 2 V[<p })- (137) 
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Therefore 

W = ^(-«) 2 + a 2 V[M)- (138) 

Now, using the definition of the slow roll parameter e = 1 — W /Ti 2 1 the expression above for H' and the second of 
Eqs.© again to give W 2 , we have 

3(^) 2 ^ 3(^) 2 



(^) 2 + a 2 U[</> ]) (a 2 ^[<M) 

where in the last step we have used the standard inflationary requirement that the potential term is much larger than 
the kinetic term. Therefore we have: (4>' ) 2 « ea 2 V/3. 

Using this result, we can rewrite the quantity c 2 appearing in Eq. (|135p as 

\2 



3 a 



GW^z'^^^Lf (MO) 
27 Mi, v Mp 



Thus, (^?) 2 » e[F/(M Pi ) 4 ]Ari. 

However as discussed in [40fl . the effect of reheating (unless some very unusual coincidences occur) is to multiply 
this result by 1/e 2 . Thus we have: 

The observations yield {^-) 2 ~ 10~ 10 . The small departure from flatness of the spectrum is taken to indicate 
e w 10 -2 . With U 1 / 4 given by the GUT scale as 10 15 GcVw 10~ 4 M P/ , and taking I 10, we conclude that 

AT is of order 10 3 >> 1. (142) 

In section [Villi following Eq. ([77|) , we used XkT = XT >> 1, so this justifies the approximation made there. 
Since we had estimated T to be of order 10 8 Mpc = 10 22 sec, we have 

A w 10" 5 MpC _1 w 10 _19 sec _1 . (143) 

Curiously, this is not far removed from the value 10 _16 sec _1 suggested by GRWflol] in their theory of instantaneous 
collapses on position, and adopted in the CSLQ theory of continuous dynamical collapse on mass-density[? ]. 

Thus for the modes of interest we have X/k in the range of 10 -2 to 10 -7 . We use this estimate to consider the 
magnitude of the sub-leading term in Eq. 



(P) 2 - (fc/2)(AT + 1 - = =), (144) 



^l + ^l + A(X/kr 



where we indeed see that it is much smaller than the leading one. 

Therefore, we have a self-consistent assignment of values for the parameters which satisfies the approximations 
made and matches the theory with the observations. 



XI. COLLAPSE ON FIELD OPERATORS 



We have presented two examples that lead to the H-Z spectrum. In one case, the collapse-generating operator 
is P ~ 7r(k) (Section IVIII[) . whose Fourier transform is the field operator 7r(x). In the other case it is X ~ y(k) 
( Section IIX[) . whose Fourier transform is the field operator y(x). In this section we shall take a look at the state 
vector evolution written in terms of these field operators. We shall give the argument in detail in the first case: the 
second case proceeds exactly similarly. 

The state vector evolution given by Eqs. ([55|) . applied to the full set of commuting collapse-generating operators we 
have discussed, is 
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= Te~ i f-T d7 >' A -f-T dr l' / + dk^[ WS (k,,,0-2A(fc)V2# s (k)] 2 -^ r dV/ + dk^[^(k,r,0-2A(fc)^* A (k)] 2 |^^ ^ 

= fe^I-T ^'"-Jxf-T *?' f+ dk[v / It« S (k,V)-2A^27fc7r s (k)] 2 --i f? T dr,' f + dk[v^A(k,V)-2A^27fc7r A (k)] 2 ,^ (3^45) 

where we have written A = A/fc and, following Eq. ([20| . we have written Ps.aO*-) = V^^s,A(k). We now define 
ws(k,t') = ^/2/kw R (k,if), w A (k 1 t') = yfzjkwi^t 1 ), and w{k,t') = w R (k,t') + iw I (k,t'). We also recall that 
7r(k) = TTs(k) + i7ivi(k). Therefore, Eq. ()145|) may be written as 

l^jy) = Te -^f^T d ^'^- 2 Jxf-T d ^'U dk l w ( k ^)~ 2 ~ Xk ' 1/2 *( k n w *^ 

= 7 - e -»/_V^'^-s/-r^Vrfk[^(k,V)-2Afc- 1 / 2 #(k)][ l0 *(k,V)-2Afc- 1 / 2 #+(k)]|^ ; _^ / 14 g^ 

where we have replaced 2x the integral by an integral which includes the lower half k-plane by defining w R (— k, t') = 
w R (k,t'), w!(-k,t r ) = -t«/(k,i / ). 

We may now convert to a real noise function and a hermitian field operator defined as 

wfrrf) = — ^ J dke^wforf), fr(x) J dke^^n(k) = (-V 2 ) -1 / 4 *(x). (147) 

Putting the inverse Fourier transforms of Eqs. (|147p into (| 146[) . we get the result we have been seeking: 

Similarly, in the second case, defining 

= J dk^k^yik) = (-V 2 ) 1 / 4 ^^, (149) 

the result is 

\1t,t) =Te- iI2 T dT >' A -kf-T d v' I rfxMx,,/)-2A5(x)] 2 |^_ r ^ ^ 15Q ^ 

We can inquire now about the fundamental reasons determining the appearance of the operators (— V 2 )~ 1 / 4 7r(x) 
(or (— V 2 ) 1 / 4 y(x)) rather than the more natural 7r(x) (or j/(x)). We believe that a truly satisfactory answer will have 
to wait for a general theory determining, in all situations, from particle physics to cosmology, the exact form of the 
CSL-type of modification to the evolution of quantum states. However we might look at the particular situation at 
hand and offer a two part response: 

First, note that the dimension of w(x',rf) and A-7r(x') have to be the same to achieve the CSL form exhibited in 
(|148[) . In order that the exponent in (|148[) be dimensionless, the dimension of 7u(x',?/) is (time) -5 / 2 , so 7r(x') must 
have the dimension (time) -3 / 2 . Now, in order that the hamiltonian in (fTU|) have the dimension (time) -1 , tt(x') has 
the dimension (time) -2 . Therefore, if the collapse-generating operator is to be of the form (operator)x7r(x'), that 
operator must have the dimension (time) 1 / 2 , which of course is the dimension of (—V 2 ) -1 / 4 . 

But, second, then it follows that the collapse-generating operator is effectively 7f(x). As we have seen, if we set 
H = in Eq. (|148|) . then at infinite time the collapse takes the system into some eigenstate of the complete commuting 
set of of operators 7r(x). But, that state will also be an eigenstate of the complete commuting set of of operators 
tt(x), since if tt(x)|0) = u(x)|0>, then tt(x)|0) = (-V 2 ) 1/4 7t(x)|0). 

Therefore, we may regard the collapse-generating operator (— \7 2 ) -1 / 4 tt(x) as the theory's way of giving us, in the 
most natural way, collapse toward eigenstates of the inflaton fluctuation momentum field 7r(x) consistent with the 
CSL state vector evolution form. It is therefore quite remarkable that such collapse leads to a prediction that agrees 
with the scale invariant spectrum, in agreement with the current cosmological observations. 

A similar case may be made for Eq. (|150l) and its effective collapse-generating operator, the inflaton fluctuation field 
y( x )- 

XII. DISCUSSION 

In this manuscript we have treated the problem of the emergence of the seeds of cosmic structure, the primordial 
inhomogeneities and anisotropics of the universe, from the quantum uncertainties (often called quantum fluctuations) 
of the inflaton field initially in the Bunch-Davies vacuum, through a version of the CSL dynamical collapse theory 
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adapted to this problem. We have taken as the basic theoretical setting the approximation known as semiclassical 
gravity. Here, gravitation is treated at the classical level while the matter fields are treated at the quantum level, and 
their energy momentum is taken to appear in Einstein's equations as the corresponding expectation value. 

In the present treatment, the expectation value of an operator differs from that given in the standard quantum 
mechanical treatment as a result of the CSL theory's modification of Schrodingcr's equation. To apply the CSL 
theory, it is necessary to select the collapse-generating operator, toward whose eigenstates the collapse occurs. We 
have considered two cases. One is where the collapse-generating operator represents the Fourier mode of the inflaton 
field perturbation. The other is where the collapse-generating operator represents the Fourier mode of the momentum 
conjugate to the inflaton field perturbation. 

As we have indicated, the quantity of direct observational interest, |a/ m | 2 , which refers to the time of the decoupling 
rjD, was evaluated instead at the end of inflation at rj = —t « 0. The changes over this interval, codified in transfer 
functions T^[r\ ri 770), when taken into account, effectively give direct observational access to the spectrum at the end 
of inflation. This turns out to be the Harrison-Zel'dovich scale-invariant spectrum. That is what the theory must 
give. 

We have seen that agreement between observations and theory will result if (P) 2 (k) turns out to be proportional 
to k. This is achieved in both cases of collapse-generating operators, by choosing the collapse rate parameter A of the 
CSL theory to have a simple dependence upon the mode's momentum magnitude k. 

In the case of operator P, it is necessary that A = A/fc. A in this case is dimensionless, so A has the dimension of 
rate. 

In the case of operator X, it is necessary that A = \k. A in this case has dimensions (time) -2 so, again, A has the 
dimension of rate. 

If we take seriously the idea that the choice of X(k) at large and small values of k must be chosen to make the 
integral in (|124|) finite, we would expect deviations in the H-Z spectral form, for large and small k that, although 
unobscrvable today, could be potentially detected in future experiments. 

Finally, one could ask which one of the two options we have discussed (or perhaps another one) might be the 
correct one? Why, in order to be consistent with observations (the H-Z spectrum), does the CSL parameter A have 
(in each case) a particular, simple, dependence on k, leading to the collapse-generating operators (— V 2 )~ 1//4 7r(x) and 
(— V 2 ) 1 / 4 y(x)? We have no deeper reason, other than "it works out that way." 

One may argue, utilizing the remarks above, that the H-Z spectrum arises from choosing A(fc) in the simplest way 
consistent with the constraints that the effective collapse parameter A ought to have the dimension of rate, that the 

evolution of each mode contains just one dimensional parameter, k, and our results for (P) 2 (k). This is hardly a 
fundamental or convincing argument. For example, in the more general context of the inflationary problem we have 
considered, there are other quantities with the same dimension as k, such as the Planck mass, the GUT scale, the 
(inflationary potential) 1 ^ 4 , the mass of the inflaton field, or a combination of these. If they are included along with 
k, there are other options to give A the dimension of rate which do not yield the H-Z spectrum. Therefore, one still 
wishes for a deeper reason for the choices we have uncovered. A response to this important question may have to wait 
for the so-far phenomenologically motivated inclusion of dynamical collapse in quantum physics to be replaced by a 
general and fundamental theory, perhaps for a natural relationship between dynamical collapse and gravitation^ . [4lj 
to be uncovered. 

We conclude that the CSL theory, with a suitable choice of the operator controlling the collapse, is capable of 
addressing the shortcomings in the standard inflationary account of the emergence of the seeds of cosmic structured 
That is, it can choose a universe possessing inhomogeneities and anisotropics. Our results are in agreement with 
observation in the regimes so far investigated empirically. It can be argued that they are also consistent with and, 
indeed, require, deviations at much smaller angular scales, which could be uncovered when the required technology 
becomes available. 



While this paper was being written, we learned of the results of a similar study that came to a different conclusion l42f . We believe that 
the reason their uncertainty, or spread of the final wave function for the relevant modes k, becomes large and ours small is connected 
with the fact that we are considering a different observable. They looked at (or took as 'focus' operator) the field amplitude for the 
v field (the Mukahnov-Sasaki variable, which is a combination of the perturbed scalar field and the Newtonian potential). We took 
as focus operator the momentum conjugate to the field variable, as discussed at the start of section VI. It is well known that the 
exponential expansion of the scale factor produces an extreme squeezing in the quantum states. That leads to an enormous growth in 
the uncertainty of the field amplitude, and an enormous decrease in the uncertainty of the conjugate momentum: the former is their 
focus operator, the latter is ours. This is consistent with our finding that the fractional dispersion in P decreases like T" 1 . In fact we 
have carried out an analysis similar to the one performed in sections VIII and IX, but taking the focus operator to be X. In that case 
wc find that the uncertainty in the value of X for the state that results from the CSL evolution diverges when the conformal time for 
the end of inflation r — > 0. 
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Appendix A: Finding w{rf) in terms of v{rj). 
In this appendix, we invert Eq. (jll3D , 

v(r)) = w(r)) - 2S J ' dT 1 'w{T 1 l )e- s{ll - 11 ' ) cos R(r] - 77'), 



(Al) 



solving for 777(77) in terms of v{rf). 

We extend the integral's lower limit to —00, keeping in mind that 777(77) = v(ri) = for 77 < —7", and extend the 
upper limit to 00 by putting a factor <d(r/ — 77') in the integrand (6(x) is the step function). Upon expressing w(rf) in 
the integrand in terms of its fourier transform w(u), and changing the integration variable to x = rf — 77, we obtain 



= w(rj) 

Taking the fourier transform yields 

v(lu) = w{ui) 



v{rf) = w(r}) -2S dx<d(~x)e Sx cosRx / due lu(x+ll) w(u) 

J —OO J — OO 

" . ium 2S(S + iw) 



{S + ilu) 2 + K 2 



or w(uj) = vluj) 5 -5 — v(oj). 

U! z — K z 
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where we have used (|109[) . When taking the fourier transform of (|A3[) . the integration path is taken below the poles 
on the real axis, so that w(rf) depends upon v(rj') for 77' < ?y. This results in 



7/7(77) = 77(7/) — 25 / drfv(r]') 

J —00 

— v ( r l) + -r f drfv(rf 
which is the result reported in Eq. (|115l) . 
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(uj — k — ie) (cj + k — ie) 
0(7/ — 7/) sin k(i] — 7/), 
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